% This file was created with JabRef 2.6.
% Encoding: UTF-8

@ARTICLE{aldunate2006,
  author = {Aldunate, Gary C AND Pestana, Reynam C},
  title = {M\~A\copyrighttodo h\~A\-brido de migra\~A\S\~A\poundso pr\~A\copyright-empilhamento
	em profundidade no dom\~A\-nio da freq\~A\textonequater\~A\textordfemeninencia
	em duas etapas com interpola\~A\S\~A\poundso},
  journal = {{Revista Brasileira de Geof\~A\-sica}},
  year = {2006},
  volume = {24},
  pages = {91 - 102},
  month = {03},
  crossref = {10.1590/S0102-261X2006000100007},
  file = {:./bibs/aldunate_pestana_2005.pdf:PDF},
  issn = {0102-261X},
  language = {pt},
  publisher = {scielo},
  url = {http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-261X2006000100007&nrm=iso}
}

@ARTICLE{mitrofanov2009,
  author = {Mitrofanov, Georgy AND Priimenko, Viatcheslav Ivanovich AND Miss\~A!`gia,
	Roseane Marchezi AND Amaral, Luis Henrique},
  title = {Transformada de Laplace na solu\~A\S\~A\poundso de problemas inversos
	din\~A\textcentmicos da s\~A\-smica},
  journal = {{Revista Brasileira de Geof\~A\-sica}},
  year = {2009},
  volume = {27},
  pages = {527 - 544},
  month = {12},
  crossref = {10.1590/S0102-261X2009000400001},
  file = {:./bibs/mitrofanov_2010.pdf:PDF},
  issn = {0102-261X},
  language = {pt},
  publisher = {scielo},
  url = {http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-261X2009000400001&nrm=iso}
}

@ARTICLE{Abubakar2008,
  author = {A Abubakar and W Hu and P M van den Berg and T M Habashy},
  title = {A finite-difference contrast source inversion method},
  journal = {Inverse Problems},
  year = {2008},
  volume = {24},
  pages = {065004},
  number = {6},
  abstract = {We present a contrast source inversion (CSI) algorithm using a finite-difference
	(FD) approach as its backbone for reconstructing the unknown material
	properties of inhomogeneous objects embedded in a known inhomogeneous
	background medium. Unlike the CSI method using the integral equation
	(IE) approach, the FD-CSI method can readily employ an arbitrary
	inhomogeneous medium as its background. The ability to use an inhomogeneous
	background medium has made this algorithm very suitable to be used
	in through-wall imaging and time-lapse inversion applications. Similar
	to the IE-CSI algorithm the unknown contrast sources and contrast
	function are updated alternately to reconstruct the unknown objects
	without requiring the solution of the full forward problem at each
	iteration step in the optimization process. The FD solver is formulated
	in the frequency domain and it is equipped with a perfectly matched
	layer (PML) absorbing boundary condition. The FD operator used in
	the FD-CSI method is only dependent on the background medium and
	the frequency of operation, thus it does not change throughout the
	inversion process. Therefore, at least for the two-dimensional (2D)
	configurations, where the size of the stiffness matrix is manageable,
	the FD stiffness matrix can be inverted using a non-iterative inversion
	matrix approach such as a Gauss elimination method for the sparse
	matrix. In this case, an LU decomposition needs to be done only once
	and can then be reused for multiple source positions and in successive
	iterations of the inversion. Numerical experiments show that this
	FD-CSI algorithm has an excellent performance for inverting inhomogeneous
	objects embedded in an inhomogeneous background medium.},
  file = {:./bibs/abubakar2008.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.06.15},
  url = {http://stacks.iop.org/0266-5611/24/i=6/a=065004}
}

@ARTICLE{Abubakar2009,
  author = {Abubakar, Aria and Hu, Wenyi and Habashy, Tarek M. and van den Berg,
	Peter M.},
  title = {Application of the finite-difference contrast-source inversion algorithm
	to seismic full-waveform data},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCC47-58},
  number = {6},
  abstract = {We have applied the finite-difference contrast-source inversion (FDCSI)
	method to seismic full-waveform inversion problems. The FDCSI method
	is an iterative nonlinear inversion algorithm. However, unlike the
	nonlinear conjugate gradient method and the Gauss-Newton method,
	FDCSI does not solve any full forward problem explicitly in each
	iterative step of the inversion process. This feature makes the method
	very efficient in solving large-scale computational problems. It
	is shown that FDCSI, with a significant lower computation cost, can
	produce inversion results comparable in quality to those produced
	by the Gauss-Newton method and better than those produced by the
	nonlinear conjugate gradient method. Another attractive feature of
	the FDCSI method is that it is capable of employing an inhomogeneous
	background medium without any extra or special effort. This feature
	is useful when dealing with time-lapse inversion problems where the
	objective is to reconstruct changes between the baseline and the
	monitor model. By using the baseline model as the background medium
	in crosswell seismic monitoring problems, high quality time-lapse
	inversion results are obtained.},
  doi = {10.1190/1.3250203},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/74/6/WCC47.pdf},
  file = {:./bibs/abubakar2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.06.15},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/74/6/WCC47}
}

@INPROCEEDINGS{Akcelik2002,
  author = {Akcelik, Volkan and Biros, George and Ghattas, Omar},
  title = {Parallel multiscale Gauss-Newton-Krylov methods for inverse wave
	propagation},
  booktitle = {Proceedings of the 2002 ACM/IEEE conference on Supercomputing},
  year = {2002},
  series = {Supercomputing '02},
  pages = {1--15},
  address = {Los Alamitos, CA, USA},
  publisher = {IEEE Computer Society Press},
  acmid = {762827},
  file = {:./bibs/akcelik2002.pdf:PDF},
  location = {Baltimore, Maryland},
  numpages = {15},
  owner = {franciane},
  timestamp = {2011.02.14},
  url = {http://portal.acm.org/citation.cfm?id=762761.762827}
}

@ARTICLE{Amestoy2000,
  author = {P. R. Amestoy and I. S. Duff and J. -Y. L'Excellent},
  title = {Multifrontal parallel distributed symmetric and unsymmetric solvers},
  journal = {Computer Methods in Applied Mechanics and Engineering},
  year = {2000},
  volume = {184},
  pages = {501 - 520},
  number = {2-4},
  abstract = {We consider the solution of both symmetric and unsymmetric systems
	of sparse linear equations. A new parallel distributed memory multifrontal
	approach is described. To handle numerical pivoting efficiently,
	a parallel asynchronous algorithm with dynamic scheduling of the
	computing tasks has been developed. We discuss some of the main algorithmic
	choices and compare both implementation issues and the performance
	of the LDLT and LU factorizations. Performance analysis on an IBM
	SP2 shows the efficiency and the potential of the method. The test
	problems used are from the Rutherford-Boeing collection and from
	the PARASOL end users.},
  doi = {DOI: 10.1016/S0045-7825(99)00242-X},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/amestoy2000.pdf:PDF},
  issn = {0045-7825},
  keywords = {MPI, MUMPS},
  owner = {franciane},
  timestamp = {2011.08.16},
  url = {http://www.sciencedirect.com/science/article/pii/S004578259900242X}
}

@ARTICLE{Amestoy2001,
  author = {Amestoy, Patrick R. and Duff, Iain S. and L'Excellent, Jean-Yves
	and Koster, Jacko},
  title = {A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic
	Scheduling},
  journal = {SIAM J. Matrix Anal. Appl.},
  year = {2001},
  volume = {23},
  pages = {15--41},
  month = {January},
  acmid = {587825},
  address = {Philadelphia, PA, USA},
  doi = {10.1137/S0895479899358194},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/amestoy2001.pdf:PDF},
  issn = {0895-4798},
  issue = {1},
  keywords = {Gaussian elimination, asynchronous parallelism, distributed memory
	computation, dynamic scheduling, multifrontal methods, sparse linear
	equations,MUMPS},
  numpages = {27},
  publisher = {Society for Industrial and Applied Mathematics},
  url = {http://portal.acm.org/citation.cfm?id=587708.587825}
}

@ARTICLE{Amestoy2006,
  author = {Patrick R. Amestoy and Abdou Guermouche and Jean-Yves L'Excellent
	and Stéphane Pralet},
  title = {Hybrid scheduling for the parallel solution of linear systems},
  journal = {Parallel Computing},
  year = {2006},
  volume = {32},
  pages = {136 - 156},
  number = {2},
  note = {Parallel Matrix Algorithms and Applications (PMAA'04)},
  abstract = {We consider the problem of designing a dynamic scheduling strategy
	that takes into account both workload and memory information in the
	context of the parallel multifrontal factorization. The originality
	of our approach is that we base our estimations (work and memory)
	on a static optimistic scenario during the analysis phase. This scenario
	is then used during the factorization phase to constrain the dynamic
	decisions that compute fully irregular partitions in order to better
	balance the workload. We show that our new scheduling algorithm significantly
	improves both the memory behaviour and the factorization time. We
	give experimental results for large challenging real-life 3D problems
	on 64 and 128 processors.},
  doi = {DOI: 10.1016/j.parco.2005.07.004},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/amestoy2006.pdf:PDF},
  issn = {0167-8191},
  keywords = {Sparse matrices, MUMPS},
  owner = {franciane},
  timestamp = {2011.08.16},
  url = {http://www.sciencedirect.com/science/article/pii/S0167819105001328}
}

@ARTICLE{Antunes2010,
  author = {Antunes, Pedro R.S. and Valtchev, Svilen S.},
  title = {A meshfree numerical method for acoustic wave propagation problems
	in planar domains with corners and cracks},
  journal = {Journal of Computational and Applied Mathematics},
  year = {2010},
  volume = {234},
  pages = {2646--2662},
  number = {9},
  month = sep,
  abstract = {The numerical solution of acoustic wave propagation problems in planar
	domains with corners and cracks is considered. Since the exact solution
	of such problems is singular in the neighborhood of the geometric
	singularities the standard meshfree methods, based on global interpolation
	by analytic functions, show low accuracy. In order to circumvent
	this issue, a meshfree modification of the method of fundamental
	solutions is developed, where the approximation basis is enriched
	by an extra span of corner adapted non-smooth shape functions. The
	high accuracy of the new method is illustrated by solving several
	boundary value problems for the Helmholtz equation, modelling physical
	phenomena from the fields of room acoustics and acoustic resonance.},
  booktitle = {Third International Workshop on Analysis and Numerical Approximation
	of Singular Problems [IWANASP08]},
  file = {:./bibs/antunes2010.pdf:PDF},
  issn = {0377-0427},
  keywords = {Meshfree methods, Method of fundamental solutions, Singular problems,
	Acoustic wave propagation, Room acoustics, Acoustic resonance},
  owner = {franciane},
  timestamp = {2011.04.04},
  url = {http://www.sciencedirect.com/science/article/B6TYH-4Y7P6VB-5/2/11cc13a112a33d54f9387a7a2633df2f}
}

@BOOK{aster2005,
  title = {Parameter Estimation and Inverse Problems},
  publisher = {Elsevier Academic Press},
  year = {2005},
  author = {Aster, Richard C. and Borchers, Brian and Thurber, Clifford H.},
  edition = {1.},
  citeulike-article-id = {4277262},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Parameter_Estimation_and_Inverse_Problems_0120656043.pdf:PDF},
  keywords = {dipl, fromarashdipl},
  posted-at = {2009-04-06 05:18:48},
  priority = {2}
}

@ARTICLE{Berenger1994,
  author = {Berenger, Jean-Pierre},
  title = {A perfectly matched layer for the absorption of electromagnetic waves},
  journal = {J. Comput. Phys.},
  year = {1994},
  volume = {114},
  pages = {185--200},
  month = {October},
  acmid = {195266},
  address = {San Diego, CA, USA},
  doi = {10.1006/jcph.1994.1159},
  file = {:./bibs/Berenger_1994.pdf:PDF},
  issn = {0021-9991},
  issue = {2},
  numpages = {16},
  publisher = {Academic Press Professional, Inc.},
  url = {http://portal.acm.org/citation.cfm?id=195261.195266}
}

@ARTICLE{vandenBerg2001,
  author = {P. M. van den Berg and A. Abubakar},
  title = {Contrast source inversion method: state of art},
  journal = {Progress In Electromagnetics Research},
  year = {2001},
  volume = {34},
  pages = {189-218},
  file = {:./bibs/vandenBerg2001.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.05.09},
  url = {http://www.jpier.org/pier/pier.php?paper=0106113}
}

@ARTICLE{Bertero2010,
  author = {M Bertero and P Boccacci and G Talenti and R Zanella and L Zanni},
  title = {A discrepancy principle for Poisson data},
  journal = {Inverse Problems},
  year = {2010},
  volume = {26},
  pages = {105004},
  number = {10},
  abstract = {In applications of imaging science, such as emission tomography, fluorescence
	microscopy and optical/infrared astronomy, image intensity is measured
	via the counting of incident particles (photons, Î³-rays, etc). Fluctuations
	in the emission-counting process can be described by modeling the
	data as realizations of Poisson random variables (Poisson data).
	A maximum-likelihood approach for image reconstruction from Poisson
	data was proposed in the mid-1980s. Since the consequent maximization
	problem is, in general, ill-conditioned, various kinds of regularizations
	were introduced in the framework of the so-called Bayesian paradigm.
	A modification of the well-known Tikhonov regularization strategy
	results in the data-fidelity function being a generalized KullbackâLeibler
	divergence. Then a relevant issue is to find rules for selecting
	a proper value of the regularization parameter. In this paper we
	propose a criterion, nicknamed discrepancy principle for Poisson
	data, that applies to both denoising and deblurring problems and
	fits quite naturally the statistical properties of the data. The
	main purpose of the paper is to establish conditions, on the data
	and the imaging matrix, ensuring that the proposed criterion does
	actually provide a unique value of the regularization parameter for
	various classes of regularization functions. A few numerical experiments
	are performed to demonstrate its effectiveness. More extensive numerical
	analysis and comparison with other proposed criteria will be the
	object of future work.},
  file = {:./bibs/bertero2010.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.06.17},
  url = {http://stacks.iop.org/0266-5611/26/i=10/a=105004}
}

@ARTICLE{biot1962,
  author = {M. A. Biot},
  title = {Mechanics of Deformation and Acoustic Propagation in Porous Media},
  journal = {Journal of Applied Physics},
  year = {1962},
  volume = {33},
  pages = {1482-1498},
  number = {4},
  doi = {10.1063/1.1728759},
  publisher = {AIP},
  url = {http://link.aip.org/link/?JAP/33/1482/1}
}

@ARTICLE{Bleibinhaus2009,
  author = {Bleibinhaus, Florian and Rondenay, Stephane},
  title = {Effects of surface scattering in full-waveform inversion},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCC69-77},
  number = {6},
  abstract = {In full-waveform inversion of seismic body waves, often the free surface
	is ignored on grounds of computational efficiency. A synthetic study
	was performed to investigate the effects of this simplification.
	In terms of size and frequency, the test model and data conform to
	a real long-offset survey of the upper crust across the San Andreas
	fault. Random fractal variations are superimposed on a background
	model with strong lateral and vertical velocity variations ranging
	from 1200 to 6800 m/s. Synthetic data were computed and inverted
	for this model and different topographies. A fully viscoelastic time-domain
	code was used to synthesize the seismograms, and a viscoacoustic
	frequency-domain code was utilized to invert them. The inversion
	was focused on early arrivals, which are dominated by P-waves but
	also contain strong P-Rayleigh wave conversions from the near-field
	of the receiver. Resulting waveform models show artifacts and a loss
	of resolution from neglecting the free surface in the inversion,
	but the inversions are stable, and they still improve the resolution
	of kinematic models. The extent of deterioration depends more on
	the subsurface than on the surface structure. Inversion results were
	improved at no additional expense by introducing a weak contrast
	along a staircase function above shots and receivers.},
  doi = {10.1190/1.3223315},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/74/6/WCC69.pdf},
  file = {:./bibs/bleibinhaus2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.14},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/74/6/WCC69}
}

@ARTICLE{,
  author = {S Bonettini and V Ruggiero},
  title = {A discrepancy principle for Poisson data: uniqueness of the solution
	for 2D and 3D data},
  year = {2010},
  file = {:./bibs/bonettini2010.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.06.17},
  url = {http://eprints.unife.it/249/}
}

@ARTICLE{Boonyasiriwat2009,
  author = {Boonyasiriwat, Chaiwoot and Valasek, Paul and Routh, Partha and Zhu,
	Xianhuai},
  title = {Application of multiscale waveform tomography for high-resolution
	velocity estimation in complex geologic environments: Canadian Foothills
	synthetic data example},
  journal = {The Leading Edge},
  year = {2009},
  volume = {28},
  pages = {454-456},
  number = {4},
  abstract = {Seismic imaging in compressional belts such as the Canadian Foothills
	is very challenging due to complex geological structures, rugged
	surface topography, and highly variable near-surface conditions.
	Seismic sections across the Canadian Foothills are usually progressively
	more distorted when approaching the Canadian Foothills region. Figure
	1 shows the degree of structural complexity and topographic variations
	which are in part responsible for the deteriorated imaging in the
	thrust belt. Accurate velocity models of subsurface structures are
	critical for improving seismic images of thrust belts in both the
	time domain (e.g., tomostatics) and the depth domain (e.g., prestack
	depth migration).},
  doi = {10.1190/1.3112764},
  eprint = {http://tle.geoscienceworld.org/cgi/reprint/28/4/454.pdf},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/Boonyasiriwat2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.14},
  url = {http://tle.geoscienceworld.org/cgi/content/abstract/28/4/454}
}

@ARTICLE{brossier2010,
  author = {Brossier, Romain and Operto, Stephane and Virieux, Jean},
  title = {Which data residual norm for robust elastic frequency-domain full
	waveform inversion?},
  journal = {Geophysics},
  year = {2010},
  volume = {75},
  pages = {R37-46},
  number = {3},
  abstract = {Elastic full-waveform inversion is an ill-posed data-fitting procedure
	that is sensitive to noise, inaccuracies of the starting model, definition
	of multiparameter classes, and inaccurate modeling of wavefield amplitudes.
	We have investigated the performance of different minimization functionals
	as the least-squares norm [IMG]f1.gif" ALT="Formula" BORDER="0">,
	the least-absolute-values norm [IMG]f2.gif" ALT="Formula" BORDER="0">,
	and combinations of both (the Huber and so-called hybrid criteria)
	with reference to two noisy offshore (Valhall model) and onshore
	(overthrust model) synthetic data sets. The four minimization functionals
	were implemented in 2D elastic frequency-domain full-waveform inversion
	(FWI), where efficient multiscale strategies were designed by successive
	inversions of a few increasing frequencies. For the offshore and
	onshore case studies, the [IMG]f2.gif" ALT="Formula" BORDER="0">-norm
	provided the most reliable models for P- and S-wave velocities ([IMG]f3.gif"
	ALT="Formula" BORDER="0"> and[IMG]f4.gif" ALT="Formula" BORDER="0">),
	even when strongly decimated data sets that correspond to few frequencies
	were used in the inversion and when outliers polluted the data. The
	[IMG]f1.gif" ALT="Formula" BORDER="0">-norm can provide reliable
	results in the presence of uniform white noise for [IMG]f3.gif" ALT="Formula"
	BORDER="0"> and [IMG]f4.gif" ALT="Formula" BORDER="0"> if the data
	redundancy is increased by refining the frequency sampling interval
	in the inversion at the expense of computational efficiency. The
	[IMG]f2.gif" ALT="Formula" BORDER="0">-norm and the Huber and hybrid
	criteria, unlike the [IMG]f1.gif" ALT="Formula" BORDER="0">-norm,
	allow for successful imaging of the [IMG]f4.gif" ALT="Formula" BORDER="0">
	model from noisy data in a soft-seabed environment, where the P-to-S-waves
	have a small footprint in the data. However, the Huber and hybrid
	criteria are sensitive to a threshold criterion that controls the
	transition between the criteria and that requires tedious trial-and-error
	investigations for reliable estimation. The [IMG]f2.gif" ALT="Formula"
	BORDER="0">-norm provides a robust alternative to the [IMG]f1.gif"
	ALT="Formula" BORDER="0">-norm for inverting decimated data sets
	in the framework of efficient frequency-domain FWI.},
  doi = {10.1190/1.3379323},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/75/3/R37.pdf},
  file = {:./bibs/brossier2010.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/75/3/R37}
}

@ARTICLE{Brossier2009,
  author = {Brossier, Romain and Operto, Stephane and Virieux, Jean},
  title = {Seismic imaging of complex onshore structures by 2D elastic frequency-domain
	full-waveform inversion},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCC105-118},
  number = {6},
  abstract = {Quantitative imaging of the elastic properties of the subsurface at
	depth is essential for civil engineering applications and oil- and
	gas-reservoir characterization. A realistic synthetic example provides
	for an assessment of the potential and limits of 2D elastic full-waveform
	inversion (FWI) of wide-aperture seismic data for recovering high-resolution
	P- and S-wave velocity models of complex onshore structures. FWI
	of land data is challenging because of the increased nonlinearity
	introduced by free-surface effects such as the propagation of surface
	waves in the heterogeneous near-surface. Moreover, the short wavelengths
	of the shear wavefield require an accurate S-wave velocity starting
	model if low frequencies are unavailable in the data. We evaluated
	different multiscale strategies with the aim of mitigating the nonlinearities.
	Massively parallel full-waveform inversion was implemented in the
	frequency domain. The numerical optimization relies on a limited-memory
	quasi-Newton algorithm thatoutperforms the more classic preconditioned
	conjugate-gradient algorithm. The forward problem is based upon a
	discontinuous Galerkin (DG) method on triangular mesh, which allows
	accurate modeling of free-surface effects. Sequential inversions
	of increasing frequencies define the most natural level of hierarchy
	in multiscale imaging. In the case of land data involving surface
	waves, the regularization introduced by hierarchical frequency inversions
	is not enough for adequate convergence of the inversion. A second
	level of hierarchy implemented with complex-valued frequencies is
	necessary and provides convergence of the inversion toward acceptable
	P- and S-wave velocity models. Among the possible strategies for
	sampling frequencies in the inversion, successive inversions of slightly
	overlapping frequency groups is the most reliable when compared to
	the more standard sequential inversion of single frequencies. This
	suggests that simultaneous inversion of multiple frequencies is critical
	when considering complex wave phenomena.},
  doi = {10.1190/1.3215771},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/74/6/WCC105.pdf},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/brossier2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.22},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/74/6/WCC105}
}

@ARTICLE{bunks1995,
  author = {Carey Bunks and Fatimetou M. Saleck and S. Zaleski and G. Chavent},
  title = {Multiscale seismic waveform inversion},
  journal = {Geophysics},
  year = {1995},
  volume = {60},
  pages = {1457-1473},
  number = {5},
  doi = {10.1190/1.1443880},
  file = {:./bibs/bunks_1995.pdf:PDF},
  publisher = {SEG},
  url = {http://link.aip.org/link/?GPY/60/1457/1}
}

@UNPUBLISHED{haroldo2008,
  author = {Haroldo Fraga de Campos Velho},
  title = {Introdução aos Problemas Inversos: Aplicações em Pesquisa Espacial},
  note = {Escola de verão em Computação Aplicada – LAC-INPE 2008},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/haroldo-Curso_PI_ELAC-2008-2.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.09.20}
}

@ARTICLE{Cerjan1985,
  author = {Charles Cerjan and Dan Kosloff and Ronnie Kosloff and Moshe Reshef},
  title = {A nonreflecting boundary condition for discrete acoustic and elastic
	wave equations},
  journal = {Geophysics},
  year = {1985},
  volume = {50},
  pages = {705-708},
  number = {4},
  doi = {10.1190/1.1441945},
  file = {:./bibs/cerjan1985.pdf:PDF},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.03.16},
  url = {http://link.aip.org/link/?GPY/50/705/1}
}

@ARTICLE{Chung2007,
  author = {Wookeen Chung and Taeyoung Ha and Wansoo Ha and Changsoo Shin},
  title = {Robust seismic waveform inversion using back-propagation algorithm},
  journal = {SEG Technical Program Expanded Abstracts},
  year = {2007},
  volume = {26},
  pages = {1780-1784},
  number = {1},
  doi = {10.1190/1.2792837},
  file = {:./bibs/chung2007.pdf:PDF},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?SGA/26/1780/1}
}

@ARTICLE{Clayton1977,
  author = {ROBERT CLAYTON and BJORN ENGQUIST},
  title = {Absorbing boundary conditions for acoustic and elastic wave equations},
  journal = {BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA},
  year = {1977},
  volume = {67},
  pages = {1529--1540},
  number = {6},
  abstract = {Boundary conditions are derived for numerical wave simulation that
	minimize artificial reflections from the edges of the domain of computation.
	In this way acoustic and elastic wave propagation in a limited area
	can be efficiently used to describe physical behavior in an unbounded
	domain. The boundary conditions are based on paraxial approximations
	of the scalar and elastic wave equations. They are computationally
	inexpensive and simple to apply, and they reduce reflections over
	a wide range of incident angles. },
  eprint = {http://www.bssaonline.org/cgi/reprint/67/6/1529.pdf},
  file = {:./bibs/Clayton1977.pdf:PDF},
  localfile = {Documents/Doutorado/bibliografia/bibs/Clayton1977.pdf},
  url = {http://www.bssaonline.org/cgi/content/abstract/67/6/1529}
}

@ARTICLE{Constable1987,
  author = {Steven C. Constable and Robert L. Parker and Catherine G. Constable},
  title = {Occam's inversion: A practical algorithm for generating smooth models
	from electromagnetic sounding data},
  journal = {Geophysics},
  year = {1987},
  volume = {52},
  pages = {289-300},
  number = {3},
  doi = {10.1190/1.1442303},
  file = {:./bibs/Constable1987.pdf:PDF},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.05.12},
  url = {http://link.aip.org/link/?GPY/52/289/1}
}

@BOOK{james1997,
  title = {Applied numerical linear algebra},
  publisher = {SIAM},
  year = {1997},
  author = {James W. Demmel},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/james_demmel_algebra.djvu:Djvu},
  owner = {franciane},
  timestamp = {2011.09.14}
}

@ARTICLE{superlu99,
  author = {James W. Demmel and Stanley C. Eisenstat and John R. Gilbert and
	Xiaoye S. Li and Joseph W. H. Liu},
  title = {A supernodal approach to sparse partial pivoting},
  journal = {SIAM J. Matrix Analysis and Applications},
  year = {1999},
  volume = {20},
  pages = {720-755},
  number = {3},
  note = {"If you use sequential SuperLU, please cite: "}
}

@ARTICLE{superlu_smp99,
  author = {James W. Demmel and John R. Gilbert and Xiaoye S. Li},
  title = {An Asynchronous Parallel Supernodal Algorithm for Sparse Gaussian
	Elimination},
  journal = {SIAM J. Matrix Analysis and Applications},
  year = {1999},
  volume = {20},
  pages = {915-952},
  number = {4},
  note = {If you use SuperLU_MT (for shared-memory parallel machines), please
	cite:}
}

@BOOK{schnabel1996,
  title = {Numerical Methods for Unconstrained Optimization and Nonlinear Equations},
  publisher = {SIAM},
  year = {1996},
  author = {J.E. Dennis and Robert B. Schnabel}
}

@BOOK{fichtner2011,
  title = {Full Seismic Waveform Modelling and Inversion},
  publisher = {Springer},
  year = {2011},
  author = {Andreas Fichtner},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Andreas Fichtner - Full seismic waveform modelling and inversion.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.06.18}
}

@BOOK{fletcher1980,
  title = {Practical Methods of Optimization - Unconstrained Op},
  publisher = {John Wiley \& Sons},
  year = {1980},
  author = {Fletcher, R.},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/fletcher.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.08.26}
}

@ARTICLE{Gallardo2007,
  author = {Gallardo, Luis A.},
  title = {Multiple cross-gradient joint inversion for geospectral imaging},
  journal = {Geophys. Res. Lett.},
  year = {2007},
  volume = {34},
  pages = {L19301--},
  number = {19},
  month = oct,
  abstract = {Accurate characterization and monitoring of complex subsurface environments
	require the integration of all the available geophysical, geochemical
	and geological information. I developed a generalized cross-gradient
	procedure that seeks multiple geometrically similar images that simplify
	the integration of cross-property subsurface information. I jointly
	invert near-surface P-wave, S-wave, DC resistivity and magnetic data
	sets recorded at a field site and compound an integrated subsurface
	(geospectral) image based on the multiple property images found.
	It is shown that, by analogy to applications in satellite imagery,
	the geospectral image assembles the multiple subsurface parameter
	values under a common structural framework that facilitates their
	visualization and analysis.},
  issn = {0094-8276},
  keywords = {joint inversion, data integration, geophysical imaging, 0902 Exploration
	Geophysics: Computational methods: seismic, 0903 Exploration Geophysics:
	Computational methods: potential fields, 0925 Exploration Geophysics:
	Magnetic and electrical methods, 3260 Mathematical Geophysics: Inverse
	theory},
  owner = {franciane},
  publisher = {AGU},
  timestamp = {2011.02.17},
  url = {http://dx.doi.org/10.1029/2007GL030409}
}

@ARTICLE{Gallardo2005,
  author = {Luis A Gallardo and Max A Meju and Marco A PÃ©rez-Flores},
  title = {A quadratic programming approach for joint image reconstruction:
	mathematical and geophysical examples},
  journal = {Inverse Problems},
  year = {2005},
  volume = {21},
  pages = {435},
  number = {2},
  abstract = {Although a comparative analysis of multiple images of a physical target
	can be useful, a joint image reconstruction approach should provide
	better interpretative elements for multi-spectral images. We present
	a generalized image reconstruction algorithm for the simultaneous
	reconstruction of band-limited images based on the novel cross-gradients
	concept developed for geophysical imaging. The general problem is
	formulated as the search for those images that stay within their
	band limits, are geometrically similar and satisfy their respective
	data in a least-squares sense. A robust iterative quadratic programming
	scheme is used to minimize the resulting objective function. We apply
	the algorithm to synthetic data generated using linear mathematical
	functions and to comparative geophysical test data. The resulting
	images recovered the test targets and show improved structural semblance
	between the reconstructed images in comparison to the results from
	two other conventional approaches.},
  file = {:./bibs/gallardo2005.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.17},
  url = {http://stacks.iop.org/0266-5611/21/i=2/a=002}
}

@ARTICLE{Gao2007,
  author = {Gao, Lingtian and Liu, Kaishin and Liu, Ying},
  title = {A meshless method for stress-wave propagation in anisotropic and
	cracked media},
  journal = {International Journal of Engineering Science},
  year = {2007},
  volume = {45},
  pages = {601--616},
  number = {2-8},
  month = feb,
  abstract = {Several transient wave propagation problems in anisotropic media and
	isotropic media with cracks are numerically analyzed by using a new
	numerical algorithm based on meshless local Petrov-Galerking (MLPG)
	method. In this algorithm, a novel modified Moving Least Squares
	(MLS) approximation is introduced to simplify the treatment of essential
	boundary conditions. By using a variant type of MLPG1 methods, the
	stabilized scheme of the discretized elasto-dynamic equations is
	obtained. Explicit central difference method with lumped mass matrix
	is used to solve the coupled ODEs to increase the efficiency of present
	algorithm. Visibility criterion is used to present the cracks, and
	the path-independent dynamic J' integral is adopted to evaluate the
	dynamic stress intensity factors. The availability and accuracy of
	the present algorithm in solving dynamic problems in isotropic or
	anisotropic media with cracks are tested through the comparison with
	the results obtained by the LS-DYNA and the method of characteristic.
	Finally, the transient stress wave interacting with a slanted crack
	under an impact loading is investigated in detail, in which the ability
	of extracting the different stress-wave components in a complex acoustic
	field is also proved.},
  issn = {0020-7225},
  keywords = {Meshless local Petrov-Galerkin method, Modified moving least square,
	Wave propagation, Dynamic stress intensity factor},
  owner = {franciane},
  timestamp = {2011.04.04},
  url = {http://www.sciencedirect.com/science/article/B6V32-4NT93Y3-1/2/e797836944fddd6e71d66a9072ac13c7}
}

@ARTICLE{Guitton2003,
  author = {Guitton, Antoine and Symes, William W.},
  title = {Robust inversion of seismic data using the Huber norm},
  journal = {Geophysics},
  year = {2003},
  volume = {68},
  pages = {1310-1319},
  number = {4},
  abstract = {The "Huber function" (or "Huber norm") is one of several robust error
	measures which interpolates between smooth ([IMG]/medium/68_4_1310_iequ1.gif"
	ALT="Formula ">2) treatment of small residuals and robust ([IMG]/medium/68_4_1310_iequ1.gif"
	ALT="Formula ">1) treatment of large residuals. Since the Huber function
	is differentiable, it may be minimized reliably with a standard gradient-based
	optimizer. We propose to minimize the Huber function with a quasi-Newton
	method that has the potential of being faster and more robust than
	conjugate-gradient methods when solving nonlinear problems. Tests
	with a linear inverse problem for velocity analysis with both synthetic
	and field data suggest that the Huber function gives far more robust
	model estimates than does a least-squares fit with or without damping.},
  doi = {10.1190/1.1598124},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/68/4/1310.pdf},
  file = {:./bibs/guitton2003.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.24},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/68/4/1310}
}

@ARTICLE{Ha2009,
  author = {Taeyoung Ha and Wookeen Chung and Changsoo Shin},
  title = {Waveform inversion using a back-propagation algorithm and a Huber
	function norm},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {R15-R24},
  number = {3},
  doi = {10.1190/1.3112572},
  file = {:./bibs/ha2009.pdf:PDF},
  keywords = {least squares approximations; oceanographic techniques; seismic waves;
	seismology; waveform analysis},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.24},
  url = {http://link.aip.org/link/?GPY/74/R15/1}
}

@ARTICLE{Abubakar2004,
  author = {T. M. Habashy and A. Abubakar},
  title = {A General Framework for Constraint Minimization for the Inversion
	of Electromagnetic Measurements},
  journal = {Progress In Electromagnetics Research},
  year = {2004},
  volume = {46},
  pages = {265-312},
  doi = {10.2528/PIER03100702},
  file = {:./bibs/Abubakar2004.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.04.27}
}

@BOOK{engl1996,
  title = {Regularization of inverse problems},
  publisher = {Kluwer Academic Publishers},
  year = {1996},
  author = {Heinz W. Engl, Martin Hanke, Andreas Neubauer},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Regularization_of_Inverse_Problems__Mathematics_and_Its_Applications_.djvu:Djvu},
  owner = {franciane},
  timestamp = {2011.09.20}
}

@ARTICLE{Herrmann2009,
  author = {Felix J. Herrmann and Yogi A. Erlangga and Tim T. Y. Lin},
  title = {Compressive simultaneous full-waveform simulation},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {A35-A40},
  number = {4},
  doi = {10.1190/1.3115122},
  file = {:./bibs/herrmann2009.pdf:PDF},
  keywords = {computational complexity; geophysics computing; seismology; waveform
	analysis},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/A35/1}
}

@ARTICLE{Hu2009,
  author = {Wenyi Hu and Aria Abubakar and Tarek M. Habashy},
  title = {Simultaneous multifrequency inversion of full-waveform seismic data},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {R1-R14},
  number = {2},
  doi = {10.1190/1.3073002},
  file = {:./bibs/hu2009.pdf:PDF},
  keywords = {finite difference time-domain analysis; geophysical techniques; Jacobian
	matrices; seismic waves; seismology},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/R1/1}
}

@ARTICLE{Hu2011,
  author = {Hu, W. and Abubakar, A. and Habashy, T. M. and Liu, J.},
  title = {Preconditioned non-linear conjugate gradient method for frequency
	domain full-waveform seismic inversion},
  journal = {Geophysical Prospecting},
  year = {2011},
  volume = {59},
  pages = {477--491},
  number = {3},
  abstract = {ABSTRACT We present preconditioned non-linear conjugate gradient algorithms
	as alternatives to the Gauss-Newton method for frequency domain full-waveform
	seismic inversion. We designed two preconditioning operators. For
	the first preconditioner, we introduce the inverse of an approximate
	sparse Hessian matrix. The approximate Hessian matrix, which is highly
	sparse, is constructed by judiciously truncating the Gauss-Newton
	Hessian matrix based on examining the auto-correlation and cross-correlation
	of the Jacobian matrix. As the second preconditioner, we employ the
	approximation of the inverse of the Gauss-Newton Hessian matrix.
	This preconditioner is constructed by terminating the iteration process
	of the conjugate gradient least-squares method, which is used for
	inverting the Hessian matrix before it converges. In our preconditioned
	non-linear conjugate gradient algorithms, the step-length along the
	search direction, which is a crucial factor for the convergence,
	is carefully chosen to maximize the reduction of the cost function
	after each iteration. The numerical simulation results show that
	by including a very limited number of non-zero elements in the approximate
	Hessian, the first preconditioned non-linear conjugate gradient algorithm
	is able to yield comparable inversion results to the Gauss-Newton
	method while maintaining the efficiency of the un-preconditioned
	non-linear conjugate gradient method. The only extra cost is the
	computation of the inverse of the approximate sparse Hessian matrix,
	which is less expensive than the computation of a forward simulation
	of one source at one frequency of operation. The second preconditioned
	non-linear conjugate gradient algorithm also significantly saves
	the computational expense in comparison with the Gauss-Newton method
	while maintaining the Gauss-Newton reconstruction quality. However,
	this second preconditioned non-linear conjugate gradient algorithm
	is more expensive than the first one.},
  doi = {10.1111/j.1365-2478.2010.00938.x},
  file = {:./bibs/hu2011.pdf:PDF},
  issn = {1365-2478},
  keywords = {Conjugate gradient, Hessian matrix, Preconditioning, Step-length},
  owner = {franciane},
  publisher = {Blackwell Publishing Ltd},
  timestamp = {2011.04.19},
  url = {http://dx.doi.org/10.1111/j.1365-2478.2010.00938.x}
}

@ARTICLE{hustedt2004,
  author = {Hustedt, Bernhard and Operto, Stéphane and Virieux, Jean},
  title = {Mixed-grid and staggered-grid finite-difference methods for frequency-domain
	acoustic wave modelling},
  journal = {Geophysical Journal International},
  year = {2004},
  volume = {157},
  pages = {1269-1296},
  number = {3},
  abstract = {SUMMARY We compare different finite-difference schemes for two-dimensional
	(2-D) acoustic frequency-domain forward modelling. The schemes are
	based on staggered-grid stencils of various accuracy and grid rotation
	strategies to discretize the derivatives of the wave equation. A
	combination of two staggered-grid stencils on the classical Cartesian
	coordinate system and the 45° rotated grid is the basis of the so-called
	mixed-grid stencil. This method is compared with a parsimonious staggered-grid
	method based on a fourth-order approximation of the first derivative
	operator. Averaging of the mass acceleration can be incorporated
	in the two stencils. Sponge-like perfectly matched layer absorbing
	boundary conditions are also examined for each stencil and shown
	to be effective. The deduced numerical stencils are examined for
	both the wavelength content and azimuthal variation. The accuracy
	of the fourth-order staggered-grid stencil is slightly superior in
	terms of phase velocity dispersion to that of the mixed-grid stencil
	when averaging of the mass acceleration term is applied to the staggered-grid
	stencil. For fourth-order derivative approximations, the classical
	staggered-grid geometry leads to a stencil that incorporates 13 grid
	nodes. The mixed-grid approach combines only nine grid nodes. In
	both cases, wavefield solutions are computed using a direct matrix
	solver based on an optimized multifrontal method. For this 2-D geometry,
	the staggered-grid strategy is significantly less efficient in terms
	of memory and CPU time requirements because of the enlarged bandwidth
	of the impedance matrix and increased number of coefficients in the
	discrete stencil. Therefore, the mixed-grid approach should be suggested
	as the routine scheme for 2-D acoustic wave propagation modelling
	in the frequency domain.},
  file = {:./bibs/hustedt_2004.pdf:PDF},
  keywords = {derivative stencils;finite-difference methods;frequency-domain wave
	modelling;seismic wave propagation}
}

@BOOK{somersalo2005,
  title = {Statistical and Computational Inverse Problems},
  publisher = {Springer},
  year = {2005},
  author = {Jari, Kaipio and Erkki, Somersalo},
  citeulike-article-id = {2396759},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/kaipio_somersallo-2005.pdf:PDF},
  posted-at = {2008-02-19 00:12:39},
  priority = {2}
}

@BOOK{nocedal,
  title = {Numerical Optimization},
  publisher = {Springer},
  year = {2006},
  author = {Jorge Nocedal, Stephen J. Wright},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Numerical_Optimization.pdf:PDF},
  keywords = {optimization},
  owner = {franciane},
  timestamp = {2011.09.05}
}

@ARTICLE{Kang2011,
  author = {Kang, Jun Won and Kallivokas, Loukas F.},
  title = {The inverse medium problem in heterogeneous PML-truncated domains
	using scalar probing waves},
  journal = {Computer Methods in Applied Mechanics and Engineering},
  year = {2011},
  volume = {200},
  pages = {265--283},
  number = {1-4},
  month = jan,
  abstract = {We discuss the inverse medium problem associated with semi-infinite
	domains. In particular, we attempt to image the spatial variability
	of shear moduli or shear wave velocities from scant surficial measurements
	of an arbitrarily heterogeneous semi-infinite domain's response to
	prescribed dynamic excitations. We use a full waveform approach to
	drive the inversion process, within a PDE-constrained optimization
	framework. Due to the semi-infinite extent of the targeted domains,
	we introduce perfectly-matched-layers (PMLs) to arrive at finite
	computational domains. The numerical implementation is based on a
	mixed finite-element method that is used to resolve the ensuing state
	and adjoint boundary-value problems, both of which are PML-endowed.
	To alleviate the inherent solution multiplicity, we use Tikhonov
	and total variation (TV) regularization schemes, in conjunction with
	a regularization factor continuation scheme. To further improve the
	optimizer's chances to converge, we also discuss a source-frequency
	continuation scheme. We report on two-dimensional numerical experiments
	using synthetic data. Included are layered profiles, and profiles
	involving inclined layers and inclusions. We also report on our methodology's
	reconstruction of the highly-heterogeneous Marmousi benchmark velocity
	model.},
  file = {:./bibs/kang2011.pdf:PDF},
  issn = {0045-7825},
  keywords = {Inverse medium problem, Full waveform inversion, Perfectly-matched-layer
	(PML), Mixed unsplit-field formulation, PDE-constrained optimization,
	Marmousi model},
  owner = {franciane},
  timestamp = {2011.04.04},
  url = {http://www.sciencedirect.com/science/article/B6V29-50X2NBP-1/2/fc282a4b60f8ba8ca72653265a2437c1}
}

@ARTICLE{Kazinnik2010,
  author = {Roman Kazinnik and Vladimir Bashkardin},
  title = {Multiscale solution of a wave equation using stable finite differences},
  journal = {SEG Technical Program Expanded Abstracts},
  year = {2010},
  volume = {29},
  pages = {3081-3086},
  number = {1},
  doi = {10.1190/1.3513487},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/kazinnik2010.pdf:PDF},
  keywords = {finite difference; migration; modeling; resolution; sampling},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?SGA/29/3081/1}
}

@ARTICLE{Komatitsch2003,
  author = {Komatitsch, Dimitri and Tromp, Jeroen},
  title = {A perfectly matched layer absorbing boundary condition for the second-order
	seismic wave equation},
  journal = {Geophysical Journal International},
  year = {2003},
  volume = {154},
  pages = {146-153},
  number = {1},
  abstract = {SUMMARY The perfectly matched layer absorbing boundary condition has
	proven to be very efficient for the elastic wave equation written
	as a first-order system in velocity and stress. We demonstrate how
	to use this condition for the same equation written as a second-order
	system in displacement. This facilitates use in the context of numerical
	schemes based upon such a system, e.g. the finite-element method,
	the spectral-element method and some finite-difference methods. We
	illustrate the efficiency of this second-order perfectly matched
	layer based upon 2-D benchmarks with body and surface waves.},
  file = {:./bibs/komatitsch_2003.pdf:PDF},
  keywords = {absorbing conditions
	
	elastic waves
	
	perfectly matched layer
	
	seismic modelling
	
	seismic wave propagation
	
	surface waves}
}

@ARTICLE{Lee2003,
  author = {Ki Ha Lee and Hee Joon Kim},
  title = {Source-independent full-waveform inversion of seismic data},
  journal = {Geophysics},
  year = {2003},
  volume = {68},
  pages = {2010-2015},
  number = {6},
  doi = {10.1190/1.1635054},
  file = {:./bibs/lee2003.pdf:PDF},
  keywords = {seismology; inverse problems; waveform analysis; frequency-domain
	analysis; geophysics computing; Earth crust},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/68/2010/1}
}

@ARTICLE{leeuwen2011,
  author = {Tristan van Leeuwen and Aleksandr Y. Aravkin and Felix J. Herrmann},
  title = {Seismic Waveform Inversion by Stochastic Optimization},
  journal = {International Journal of Geophysics},
  year = {2011},
  volume = {2011},
  pages = {18},
  doi = {doi:10.1155/2011/689041},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/leeuwen2011.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.07.04}
}

@ARTICLE{Leeuwen2010,
  author = {T van Leeuwen and W A Mulder},
  title = {A comparison of seismic velocity inversion methods for layered acoustics},
  journal = {Inverse Problems},
  year = {2010},
  volume = {26},
  pages = {015008},
  number = {1},
  abstract = {In seismic imaging, one tries to infer the medium properties of the
	subsurface from seismic reflection data. These data are the result
	of an active source experiment, where an explosive source and an
	array of receivers are placed at the surface. Because of the absence
	of low frequencies in the data, the corresponding inverse problem
	is strongly nonlinear in the slowly varying component of the velocity.
	The least-squares misfit functional typically exhibits local minima
	and has a small basin of attraction. The usual approach to fitting
	the data in a least-squares sense by employing a gradient-based optimization
	method will therefore most likely result in a wrong velocity model.
	In the geophysical community, this problem has long been recognized
	and alternative formulations of the inverse problem have been developed.
	We review several of these formulations and analyse the sensitivity
	to the error in the smooth velocity component. This analysis is carried
	out for laterally homogeneous velocities using an asymptotic solution
	of the wave equation. The analysis suggests that formulations which
	are geared towards fitting the phases of the data, rather than the
	amplitudes, have smooth corresponding misfit functionals with a large
	basin of attraction.},
  file = {:./bibs/leeuwen2010.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.14},
  url = {http://stacks.iop.org/0266-5611/26/i=1/a=015008}
}

@ARTICLE{demmel2003,
  author = {Xiaoye S. Li and James W. Demmel},
  title = {{SuperLU_DIST}: A Scalable Distributed-Memory Sparse Direct Solver
	for Unsymmetric Linear Systems},
  journal = {ACM Trans. Mathematical Software},
  year = {2003},
  volume = {29},
  pages = {110-140},
  number = {2},
  month = {June}
}

@ARTICLE{LLi2009,
  author = {Lianlin Li, Hu Zheng, and Fang Li},
  title = {Two-Dimensional Contrast Source Inversion Method With Phaseless Data:
	TM Case},
  journal = {IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING},
  year = {2009},
  volume = {47},
  pages = {1719-1736},
  abstract = {In this paper, two new approaches are presented for the solution of
	electromagnetic inverse scattering problems when amplitude-only data
	are available. The proposed techniques are based on a customized
	version, which are the so-called contrast source inversion (CSI)
	and multiplicative regularized CSI (MRCSI) methods. In the proposed
	approaches, denoted as the phaseless-data (PD)-CSI and the PD-MRCSI,
	only the term of the cost functional concerning the mismatch between
	the measured and estimated field data (i.e., the data equation) has
	been properly redefined. Moreover, the back-projection algorithm
	has been modified to provide an initial solution ensuring the rapid
	convergence of the optimization procedures and avoid the reconstruction
	of false solutions. A set of representative results concerning numerical
	as well as experimental tests is reported to show the accuracy of
	the proposed amplitude-only reconstruction approaches.},
  doi = {10.1109/TGRS.2008.2006360},
  file = {:./bibs/LLi2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.05.09},
  url = {http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=4695996&isnumber=4939375}
}

@ARTICLE{Liao1996,
  author = {Liao, Qingbo and McMechan, George A.},
  title = {Multifrequency viscoacoustic modeling and inversion},
  journal = {Geophysics},
  year = {1996},
  volume = {61},
  pages = {1371-1378},
  number = {5},
  abstract = {Modeling and inversion for seismic wavefields that include the attenuation
	and phase dispersion effects of Q can be implemented in the space-frequency
	domain. The viscoacoustic wave equation is solved by the moment method.
	Absorbing boundary conditions are implemented by reducing Q and adjusting
	the complex velocity (to reduce Q-dependent reflectivity) in a zone
	around the edges of the model grid. Nonlinear inversion is performed
	using iterative linearized inversions. The residual wavefield at
	a single frequency is back projected, using an anticausal Green's
	function, along the viscoacoustic wavepath in an estimate of the
	model, to get updated velocity and Q distributions. The model obtained
	from data at one frequency becomes input to inversion at the next
	higher frequency. Velocity and Q are inverted simultaneously as they
	are interdependent. Both modeling and inversion algorithms are successfully
	tested with synthetic examples; data at two or three frequencies
	are sufficient to produce reliable images from noise-free synthetic
	data.},
  doi = {10.1190/1.1444060},
  file = {:./bibs/LiaoMcMechan1996.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.15},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/61/5/1371}
}

@ARTICLE{Liu1989,
  author = {Liu, Dong C. and Nocedal, Jorge},
  title = {On the limited memory BFGS method for large scale optimization},
  journal = {Mathematical Programming},
  year = {1989},
  volume = {45},
  pages = {503-528},
  note = {10.1007/BF01589116},
  abstract = {We study the numerical performance of a limited memory quasi-Newton
	method for large scale optimization, which we call the L-BFGS method.
	We compare its performance with that of the method developed by Buckley
	and LeNir (1985), which combines cycles of BFGS steps and conjugate
	direction steps. Our numerical tests indicate that the L-BFGS method
	is faster than the method of Buckley and LeNir, and is better able
	to use additional storage to accelerate convergence. We show that
	the L-BFGS method can be greatly accelerated by means of a simple
	scaling. We then compare the L-BFGS method with the partitioned quasi-Newton
	method of Griewank and Toint (1982a). The results show that, for
	some problems, the partitioned quasi-Newton method is clearly superior
	to the L-BFGS method. However we find that for other problems the
	L-BFGS method is very competitive due to its low iteration cost.
	We also study the convergence properties of the L-BFGS method, and
	prove global convergence on uniformly convex problems.},
  affiliation = {Department of Electrical Engineering and Computer Science Northwestern
	University 60208 Evanston IL USA},
  issn = {0025-5610},
  issue = {1},
  keyword = {Mathematics and Statistics},
  owner = {franciane},
  publisher = {Springer Berlin / Heidelberg},
  timestamp = {2011.09.05},
  url = {http://dx.doi.org/10.1007/BF01589116}
}

@ARTICLE{HuoLiu1997,
  author = {Qing-Huo Liu and Jianping Tao},
  title = {The perfectly matched layer for acoustic waves in absorptive media},
  journal = {The Journal of the Acoustical Society of America},
  year = {1997},
  volume = {102},
  pages = {2072-2082},
  number = {4},
  doi = {10.1121/1.419657},
  file = {:./bibs/HuoLiu_1997.pdf:PDF},
  keywords = {acoustic wave propagation; acoustic wave absorption; finite difference
	time-domain analysis; underwater sound},
  publisher = {ASA},
  url = {http://link.aip.org/link/?JAS/102/2072/1}
}

@ARTICLE{Ma2010,
  author = {Yong Ma and Dave Hale and Zhaobo (Joe) Meng and Bin Gong},
  title = {Full waveform inversion with image-guided gradient},
  journal = {SEG Technical Program Expanded Abstracts},
  year = {2010},
  volume = {29},
  pages = {1003-1007},
  number = {1},
  doi = {10.1190/1.3513016},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/ma2010.pdf:PDF},
  keywords = {imaging; interpolation; inversion; velocity},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?SGA/29/1003/1}
}

@ARTICLE{Malinowski2011,
  author = {Malinowski, M. and Operto, S. and Ribodetti, A.},
  title = {High-resolution seismic attenuation imaging from wide-aperture onshore
	data by visco-acoustic frequency-domain full-waveform inversion},
  journal = {Geophysical Journal International},
  year = {2011},
  pages = {no--no},
  abstract = {SUMMARY Here we assess the potential of the visco-acoustic frequency
	domain full-waveform inversion (FWI) to reconstructâP-wave velocity
	(VP) andâP-wave attenuation factor (Q) from surface onshore seismic
	data. First, we perform a sensitivity analysis of the FWI based upon
	a grid search analysis of the misfit function and several synthetic
	FWI examples using velocity andâQâmodels of increasing complexity.
	Subsequently, we applied both the acoustic and visco-acoustic FWI
	to real surface wide-aperture onshore seismic data from the Polish
	Basin, where a strong attenuation of the seismic data is observed.
	The sensitivity analysis of the visco-acoustic FWI suggests that
	the FWI can jointly reconstruct the velocity and the attenuation
	factor if the signature of the attenuation is sufficiently strong
	in the data. A synthetic example corresponding to a homogeneous background
	model with an inclusion shows a reliable reconstruction ofâVPâandâQâin
	the inclusion, whenâQâis as small as 90 and 50 in the background
	model and in the inclusion, respectively. A first application of
	acoustic FWI to real data shows that a heuristic normalization of
	the data with offset allows us to compensate for the effect of the
	attenuation in the data and reconstruct a reliable velocity model.
	Alternatively, we show that visco-acoustic FWI allows us to reconstruct
	jointly both a reliable velocity model and aâQâmodel from the
	true-amplitude data. We propose a pragmatical approach based upon
	seismic modelling and source wavelet estimation to infer the best
	starting homogeneousâQâmodel for visco-acoustic FWI. We find
	the source wavelet estimation quite sensitive to the quality of the
	velocity and attenuation models used for the estimation. For example,
	source-to-source wavelets are significantly more consistent when
	computed in the final FWI model than in the initial one. A good kinematic
	and amplitude match between the early-arriving phases of the real
	and time-domain synthetic seismograms computed in the final FWI model
	provides an additional evidence of the reliability of the final FWI
	model. We find the recovered velocity and attenuation models consistent
	with the expected lithology and stratigraphy in the study area. We
	link high-attenuation zones with the increased clay content and the
	presence of the mineralized fluids.},
  doi = {10.1111/j.1365-246X.2011.05098.x},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/Malinowski2011.pdf:PDF},
  issn = {1365-246X},
  keywords = {Numerical solutions, Inverse theory, Controlled source seismology,
	Seismic attenuation, Seismic tomography},
  owner = {franciane},
  publisher = {Blackwell Publishing Ltd},
  timestamp = {2011.07.21},
  url = {http://dx.doi.org/10.1111/j.1365-246X.2011.05098.x}
}

@ARTICLE{martin_komatitsch_2008,
  author = {Roland Martin and Dimitri Komatitsch and Abdelaziz Ezziani},
  title = {An unsplit convolutional perfectly matched layer improved at grazing
	incidence for seismic wave propagation in poroelastic media},
  journal = {Geophysics},
  year = {2008},
  volume = {73},
  pages = {T51-T61},
  number = {4},
  doi = {10.1190/1.2939484},
  file = {:./bibs/martin_komatitsch_2008.pdf:PDF},
  keywords = {elastic waves; finite difference methods; geophysical techniques;
	wave equations},
  publisher = {SEG},
  url = {http://link.aip.org/link/?GPY/73/T51/1}
}

@BOOK{menke,
  title = {Geophysical data analysis: Discrete inverse theory},
  publisher = {Academic Press, Inc.},
  year = {1989},
  author = {William Menke},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Geophysical_Data_Analysis.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.09.05}
}

@ARTICLE{Mulder2004,
  author = {Mulder, W. A. and Plessix, R.-E.},
  title = {How to choose a subset of frequencies in frequency-domain finite-difference
	migration},
  journal = {Geophysical Journal International},
  year = {2004},
  volume = {158},
  pages = {801--812},
  number = {3},
  abstract = {SUMMARY Finite-difference migration with the two-way wave equation
	can be accelerated by an order of magnitude if the frequency domain
	rather than the time domain is used. This gain is mainly accomplished
	by using a subset of the available frequencies. The implicit assumption
	is that the data have a certain amount of redundancy in the frequency
	domain.The choice of frequencies cannot be arbitrary. If the frequencies
	are chosen with a constant increment and their spacing is too large,
	the well-known wrap-around that occurs when transforming back to
	the time domain will also show up in the migration to the depth domain,
	albeit in a more subtle way. Because migration involves propagation
	in a given background velocity model and summation over shots and
	receivers, the effects of wrap-around may disappear even when the
	Nyquist theorem is not obeyed.We have studied these effects analytically
	for the constant-velocity case and determined sampling conditions
	that avoid wrap-around artefacts. The conditions depend on the velocity,
	depth of the migration grid and offset range. They show that the
	spacing between subsequent frequencies can be larger than the inverse
	of the time range prescribed by the Nyquist theorem. A 2-D example
	has been used to test the validity of these conditions for a more
	realistic velocity model. Finite-difference migration with the one-way
	wave equation shows a similar behaviour.},
  doi = {10.1111/j.1365-246X.2004.02336.x},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/mulder2004.pdf:PDF},
  issn = {1365-246X},
  keywords = {finite-differences, frequency domain, migration},
  owner = {franciane},
  publisher = {Blackwell Science Ltd},
  timestamp = {2011.02.14},
  url = {http://dx.doi.org/10.1111/j.1365-246X.2004.02336.x}
}

@ARTICLE{Neklyudov2010,
  author = {Dmitry Neklyudov and Ilya Silvestrov and Vladimir Tcheverda},
  title = {A Helmholtz iterative solver with semianalytical preconditioner for
	the frequency-domain full-waveform inversion},
  journal = {SEG Technical Program Expanded Abstracts},
  year = {2010},
  volume = {29},
  pages = {1070-1074},
  number = {1},
  doi = {10.1190/1.3513031},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/Neklyudov2010.pdf:PDF},
  keywords = {inversion; nonlinear; VSP},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?SGA/29/1070/1}
}

@ARTICLE{nocedal1980,
  author = {Nocedal, Jorge},
  title = {Updating Quasi-Newton Matrices with Limited Storage},
  journal = {Mathematics of Computation},
  year = {1980},
  volume = {35},
  pages = {pp. 773-782},
  number = {151},
  abstract = {We study how to use the BFGS quasi-Newton matrices to precondition
	minimization methods for problems where the storage is critical.
	We give an update formula which generates matrices using information
	from the last $m$ iterations, where $m$ is any number supplied by
	the user. The quasi-Newton matrix is updated at every iteration by
	dropping the oldest information and replacing it by the newest information.
	It is shown that the matrices generated have some desirable properties.
	The resulting algorithms are tested numerically and compared with
	several well-known methods.},
  copyright = {Copyright © 1980 American Mathematical Society},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/nocedal1980.pdf:PDF},
  issn = {00255718},
  jstor_articletype = {research-article},
  jstor_formatteddate = {Jul., 1980},
  language = {English},
  publisher = {American Mathematical Society},
  url = {http://www.jstor.org/stable/2006193}
}

@ARTICLE{operto_2007,
  author = {Stephane Operto and Jean Virieux and Patrick Amestoy and Jean-Yves
	L'Excellent and Luc Giraud and Hafedh Ben Hadj Ali},
  title = {3D finite-difference frequency-domain modeling of visco-acoustic
	wave propagation using a massively parallel direct solver: A feasibility
	study},
  journal = {Geophysics},
  year = {2007},
  volume = {72},
  pages = {SM195-SM211},
  number = {5},
  doi = {10.1190/1.2759835},
  file = {:./bibs/operto_2007.pdf:PDF},
  keywords = {acoustic wave propagation; finite difference methods; geophysical
	techniques; seismic waves; seismology; wave equations},
  publisher = {SEG},
  url = {http://link.aip.org/link/?GPY/72/SM195/1}
}

@ARTICLE{operto_2009,
  author = {Stephane Operto and Jean Virieux and A. Ribodetti and J. E. Anderson},
  title = {Finite-difference frequency-domain modeling of viscoacoustic wave
	propagation in 2D tilted transversely isotropic (TTI) media},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {T75-T95},
  number = {5},
  doi = {10.1190/1.3157243},
  file = {:./bibs/operto_2009.pdf:PDF},
  keywords = {anisotropic media; seismic waves; viscoelasticity; wave equations;
	wave propagation},
  publisher = {SEG},
  url = {http://link.aip.org/link/?GPY/74/T75/1}
}

@ARTICLE{Plessix2009,
  author = {Rene-Edouard Plessix},
  title = {Three-dimensional frequency-domain full-waveform inversion with an
	iterative solver},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCC149-WCC157},
  number = {6},
  doi = {10.1190/1.3211198},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/plessix2009.pdf:PDF},
  keywords = {data acquisition; seafloor phenomena; seismic waves; seismology},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/WCC149/1}
}

@ARTICLE{plessix2006,
  author = {Plessix, R.-E.},
  title = {A review of the adjoint-state method for computing the gradient of
	a functional with geophysical applications},
  journal = {Geophysical Journal International},
  year = {2006},
  volume = {167},
  pages = {495--503},
  number = {2},
  abstract = {Estimating the model parameters from measured data generally consists
	of minimizing an error functional. A classic technique to solve a
	minimization problem is to successively determine the minimum of
	a series of linearized problems. This formulation requires the Fréchet
	derivatives (the Jacobian matrix), which can be expensive to compute.
	If the minimization is viewed as a non-linear optimization problem,
	only the gradient of the error functional is needed. This gradient
	can be computed without the Fréchet derivatives. In the 1970s, the
	adjoint-state method was developed to efficiently compute the gradient.
	It is now a well-known method in the numerical community for computing
	the gradient of a functional with respect to the model parameters
	when this functional depends on those model parameters through state
	variables, which are solutions of the forward problem. However, this
	method is less well understood in the geophysical community. The
	goal of this paper is to review the adjoint-state method. The idea
	is to define some adjoint-state variables that are solutions of a
	linear system. The adjoint-state variables are independent of the
	model parameter perturbations and in a way gather the perturbations
	with respect to the state variables. The adjoint-state method is
	efficient because only one extra linear system needs to be solved.Several
	applications are presented. When applied to the computation of the
	derivatives of the ray trajectories, the link with the propagator
	of the perturbed ray equation is established.},
  doi = {10.1111/j.1365-246X.2006.02978.x},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/plessix2006.pdf:PDF},
  issn = {1365-246X},
  keywords = {adjoint state, gradient, migration, tomography},
  publisher = {Blackwell Publishing Ltd},
  url = {http://dx.doi.org/10.1111/j.1365-246X.2006.02978.x}
}

@ARTICLE{pratt1999a,
  author = {R. Gerhard Pratt},
  title = {Seismic waveform inversion in the frequency domain, Part 1: Theory
	and verification in a physical scale model},
  journal = {Geophysics},
  year = {1999},
  volume = {64},
  pages = {888--901},
  number = {3},
  doi = {10.1190/1.1444597},
  file = {:./bibs/pratt1999.pdf:PDF},
  publisher = {SEG},
  url = {http://link.aip.org/link/?GPY/64/888/1}
}

@ARTICLE{pratt1990,
  author = {Pratt, R. Gerhard},
  title = {Frequency-domain elastic wave modeling by finite differences; a tool
	for crosshole seismic imaging},
  journal = {Geophysics},
  year = {1990},
  volume = {55},
  pages = {626-632},
  number = {5},
  abstract = {The migration, imaging, or inversion of wide-aperture cross-hole data
	depends on the ability to model wave propagation in complex media
	for multiple source positions. Computational costs can be considerably
	reduced in frequency-domain imaging by modeling the frequency-domain
	steady-state equations, rather than the time-domain equations of
	motion. I develop a frequency-domain approach in this note that is
	competitive with time-domain modeling when solutions for multiple
	sources are required or when only a limited number of frequency components
	of the solution are required.},
  doi = {10.1190/1.1442874},
  file = {:./bibs/pratt_modelling1990.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/55/5/626}
}

@ARTICLE{pratt1999b,
  author = {Pratt, R. Gerhard and Shipp, Richard M.},
  title = {Seismic waveform inversion in the frequency domain; Part 2; Fault
	delineation in sediments using crosshole data},
  journal = {Geophysics},
  year = {1999},
  volume = {64},
  pages = {902-914},
  number = {3},
  abstract = {A crosshole experiment was carried out in a layered sedimentary environment
	in which a normal fault is known to cut through the section. Initial
	traveltime inversions produced stable but low-resolution images from
	which the fault could be only vaguely inferred. To image the fault,
	wavefield inversion was used to produce a velocity model consistent
	with the detailed phase and amplitude of the data at a number of
	frequencies. Our wavefield inversion scheme uses a classical, descent-type
	algorithm for decreasing the data misfit by iteratively computing
	the gradient of this misfit by repeated forward and backward propagations.
	Our propagator is a full-wave equation, frequency-domain, acoustic,
	finite-difference method. The use of the frequency-space domain yields
	computational advantages for multisource data and allows an easy
	incorporation of viscous effects. By running wavefield inversion
	on the field data, a quantitative velocity image was produced that
	yielded a significantly improved image of the fault (when compared
	with the original traveltime inversions). Because the original field
	data were noisy and contained a high degree of multiple scattering
	(from the layering of the sediments), the transmitted arrivals were
	selectively windowed to enhance the image. The sediments at the site
	were strongly attenuating; we therefore used a viscoacoustic model
	during the modeling and the inversion that correctly simulated the
	observed decrease in amplitude with increasing frequency and source-receiver
	offset. Furthermore, since the traveltime inversion indicated a high
	degree of anisotropy at the site, a fixed, homogeneous level of anisotropy
	was used during the inversion. Tests at varying levels of anisotropy
	confirmed the improvement in image quality and in data fit when anisotropy
	was incorporated. The final image was verified by examining the distribution
	of the residuals in the frequency domain, by comparing time-domain
	modeled wavefields with the observed data, and by direct comparison
	with borehole logs.},
  doi = {10.1190/1.1444598},
  file = {:./bibs/prattshipp1999.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.14},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/64/3/902}
}

@BOOK{Ramm,
  title = {Inverse problems - Mathematical and analytical techniques with applications
	to engineering.},
  publisher = {Springer},
  year = {2005},
  author = {Alexander G. Ramm},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Inverse Problems Mathematical and Analytical Techniques with Applications to Engineering_ISBN0387231951.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.08.26}
}

@ARTICLE{Ravaut2004,
  author = {Ravaut, C. and Operto, S. and Improta, L. and Virieux, J. and Herrero,
	A. and Dell'Aversana, P.},
  title = {Multiscale imaging of complex structures from multifold wide-aperture
	seismic data by frequency-domain full-waveform tomography: application
	to a thrust belt},
  journal = {Geophysical Journal International},
  year = {2004},
  volume = {159},
  pages = {1032-1056},
  number = {3},
  abstract = {SUMMARY An application of full-waveform tomography to dense onshore
	wide-aperture seismic data recorded in a complex geological setting
	(thrust belt) is presented. The waveform modelling and tomography
	are implemented in the frequency domain. The modelling part is solved
	with a finite-difference method applied to the visco-acoustic wave
	equation. The inversion is based on a local gradient method. Only
	the P-wave velocity is involved in the inversion. The inversion is
	applied iteratively to discrete frequency components by proceeding
	from low to high frequencies. This defines a multiscale imaging in
	the sense that high wavenumbers are progressively incorporated in
	images. The linearized waveform tomography requires an accurate starting
	velocity model that has been developed by first-arrival traveltime
	tomography. After specific pre-processing of the data, 16 frequency
	components ranging between 5.4 and 20 Hz were inverted. Ten iterations
	were computed per frequency component leading to 160 tomographic
	models. The waveform tomography has successfully imaged southwest-dipping
	structures previously identified from other geophysical data as being
	associated with high-resistivity bodies. The relevance of the tomographic
	images is locally demonstrated by comparison of a velocity–depth
	function extracted from the waveform tomography models with a coincident
	vertical seismic profiling (VSP) log available on the profile. Moreover,
	comparison between observed and synthetic seismograms computed in
	the (starting) traveltime and waveform tomography models demonstrates
	unambiguously that the waveform tomography successfully predicts
	for wide-angle reflections from southwest-dipping geological structures.
	This study demonstrates that the combination of first-arrival traveltime
	and frequency-domain full-waveform tomographies applied to dense
	wide-aperture seismic data is a promising approach to quantitative
	imaging of complex geological structures. Indeed, wide-aperture acquisition
	geometries offer the opportunity to develop an accurate background
	velocity model for the subsequent waveform tomography. This is critical,
	because the building of the macromodel remains an open question when
	only near-vertical reflection data are considered.},
  file = {:./bibs/ravaut_2004.pdf:PDF},
  keywords = {finite-difference methods
	
	thrust belt
	
	traveltime and full waveform inversions
	
	wide-aperture seismic data}
}

@ARTICLE{Rudin1992,
  author = {Rudin, Leonid I. and Osher, Stanley and Fatemi, Emad},
  title = {Nonlinear total variation based noise removal algorithms},
  journal = {Phys. D},
  year = {1992},
  volume = {60},
  pages = {259--268},
  month = {November},
  acmid = {142312},
  address = {Amsterdam, The Netherlands, The Netherlands},
  doi = {http://dx.doi.org/10.1016/0167-2789(92)90242-F},
  file = {:./bibs/Rudin1992.pdf:PDF},
  issn = {0167-2789},
  issue = {1-4},
  numpages = {10},
  publisher = {Elsevier Science Publishers B. V.},
  url = {http://dx.doi.org/10.1016/0167-2789(92)90242-F}
}

@ARTICLE{Ryan1994,
  author = {Harold Ryan},
  title = {Ricker, Ormsby, Klauder, Butterworth - a choice of wavelets},
  journal = {CSEG Recorder},
  year = {1994},
  file = {:./bibs/ryan1994.pdf:PDF;:./bibs/choice-of-wavelets.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.04.04}
}

@BOOK{sheriff,
  title = {Encyclopedic Dictionary of Applied Geophysics},
  publisher = {SEG},
  year = {2002},
  author = {Robert Sheriff},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Sheriff-dic.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.08.26}
}

@ARTICLE{Shin1995,
  author = {Changsoo Shin},
  title = {Sponge boundary condition for frequency-domain modeling},
  journal = {Geophysics},
  year = {1995},
  volume = {60},
  pages = {1870-1874},
  number = {6},
  doi = {10.1190/1.1443918},
  file = {:./bibs/shin1995.pdf:PDF},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.03.16},
  url = {http://link.aip.org/link/?GPY/60/1870/1}
}

@ARTICLE{shin2001,
  author = {Shin, Changsoo and Jang, Seonghyung and Min, Dong-Joo},
  title = {Improved amplitude preservation for prestack depth migration by inverse
	scattering theory},
  journal = {Geophysical Prospecting},
  year = {2001},
  volume = {49},
  pages = {592--606},
  number = {5},
  abstract = {A prestack reverse time-migration image is not properly scaled with
	increasing depth. The main reason for the image being unscaled is
	the geometric spreading of the wavefield arising during the back-propagation
	of the measured data and the generation of the forward-modelled wavefields.
	This unscaled image can be enhanced by multiplying the inverse of
	the approximate Hessian appearing in the Gauss–Newton optimization
	technique. However, since the approximate Hessian is usually too
	expensive to compute for the general geological model, it can be
	used only for the simple background velocity model.We show that the
	pseudo-Hessian matrix can be used as a substitute for the approximate
	Hessian to enhance the faint images appearing at a later time in
	the 2D prestack reverse time-migration sections. We can construct
	the pseudo-Hessian matrix using the forward-modelled wavefields (which
	are used as virtual sources in the reverse time migration), by exploiting
	the uncorrelated structure of the forward-modelled wavefields and
	the impulse response function for the estimated diagonal of the approximate
	Hessian. Although it is also impossible to calculate directly the
	inverse of the pseudo-Hessian, when using the reciprocal of the pseudo-Hessian
	we can easily obtain the inverse of the pseudo-Hessian. As examples
	supporting our assertion, we present the results obtained by applying
	our method to 2D synthetic and real data collected on the Korean
	continental shelf.},
  doi = {10.1046/j.1365-2478.2001.00279.x},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/shin2001.pdf:PDF},
  issn = {1365-2478},
  publisher = {Blackwell Science Ltd},
  url = {http://dx.doi.org/10.1046/j.1365-2478.2001.00279.x}
}

@ARTICLE{Shin2010,
  author = {Shin, Changsoo and Koo, Nam-Hyung and Cha, Young Ho and Park, Keun-Pil
	and Ravaut, C. and Operto, S. and Improta, L. and Virieux, J. and
	Herrero, A. and Dell'Aversana, P.},
  title = {Sequentially ordered single-frequency 2-D acoustic waveform inversion
	in the Laplace–Fourier domain
	
	Multiscale imaging of complex structures from multifold wide-aperture
	seismic data by frequency-domain full-waveform tomography: application
	to a thrust belt},
  journal = {Geophysical Journal International},
  year = {2010},
  volume = {181},
  pages = {935-950},
  number = {2},
  abstract = {SUMMARY In the conventional frequency-domain waveform inversion, either
	multifrequency simultaneous inversion or sequential single-frequency
	inversion has been implemented. However, most conventional frequency-domain
	waveform inversion methods fail to recover background velocity when
	low-frequency information is missing. Recently, new waveform inversion
	techniques in the Laplace and Laplace–Fourier domain have been proposed
	to recover background velocity structure from data with insufficient
	low-frequency information. In such techniques, however, all frequencies
	are inverted simultaneously, and this requires large computational
	resources and long computation times. In this paper, we propose a
	sequentially ordered single-frequency 2-D acoustic waveform inversion
	using the logarithmic objective function in the Laplace–Fourier domain.
	Our algorithm sequentially inverts single-frequency data in the Laplace–Fourier
	domain, thus reducing computational resources. Unlike most conventional
	waveform inversion methods requiring an initial velocity model close
	to the true model, we propose a one-step waveform inversion method
	in seeking to find a final velocity structure from the simple initial
	model through a hybrid combination of the Laplace domain inversion
	and the Fourier domain inversion. We adopt and evaluate the multiloop
	algorithm by modifying the double-loop algorithm commonly used in
	the conventional frequency-domain waveform inversion. Using the multiloop
	algorithm repeating loop over frequencies, the quality of the inversion
	results can be improved and the decision problem of the number of
	iterations for each frequency can be overcome effectively. Because
	the sequential order of the Laplace–Fourier frequencies in a 2-D
	plane should be assigned for inverting Laplace–Fourier frequency
	data consecutively, we present three different sequential orders
	of Laplace–Fourier frequencies while considering the multiscale and
	layer-stripping approach, and we compare the inversion results from
	the numerical experiments. We applied the sequentially ordered single-frequency
	2-D acoustic waveform inversion in the full Laplace–Fourier domain
	to the synthetic seismic data produced from complex structure model
	and field data. A realistic model could be recovered in an efficient
	and robust manner, even using the two-layer homogeneous velocity
	model as an initial model. The inverted velocity model from the field
	data was validated by examining the migrated image from the pre-stack
	depth migration and the flattening of the common-image gathers or
	by comparing the synthetic shot gather with the real shot gather.
	The proposed one-step waveform inversion algorithm can be easily
	extended to the sequential inversion of 3-D acoustic or elastic data
	in the full Laplace–Fourier domain. ER Provider: John Wiley & Sons,
	Ltd Content:text/plain; charset="UTF-8"
	
	SUMMARY An application of full-waveform tomography to dense onshore
	wide-aperture seismic data recorded in a complex geological setting
	(thrust belt) is presented. The waveform modelling and tomography
	are implemented in the frequency domain. The modelling part is solved
	with a finite-difference method applied to the visco-acoustic wave
	equation. The inversion is based on a local gradient method. Only
	the P-wave velocity is involved in the inversion. The inversion is
	applied iteratively to discrete frequency components by proceeding
	from low to high frequencies. This defines a multiscale imaging in
	the sense that high wavenumbers are progressively incorporated in
	images. The linearized waveform tomography requires an accurate starting
	velocity model that has been developed by first-arrival traveltime
	tomography. After specific pre-processing of the data, 16 frequency
	components ranging between 5.4 and 20 Hz were inverted. Ten iterations
	were computed per frequency component leading to 160 tomographic
	models. The waveform tomography has successfully imaged southwest-dipping
	structures previously identified from other geophysical data as being
	associated with high-resistivity bodies. The relevance of the tomographic
	images is locally demonstrated by comparison of a velocity–depth
	function extracted from the waveform tomography models with a coincident
	vertical seismic profiling (VSP) log available on the profile. Moreover,
	comparison between observed and synthetic seismograms computed in
	the (starting) traveltime and waveform tomography models demonstrates
	unambiguously that the waveform tomography successfully predicts
	for wide-angle reflections from southwest-dipping geological structures.
	This study demonstrates that the combination of first-arrival traveltime
	and frequency-domain full-waveform tomographies applied to dense
	wide-aperture seismic data is a promising approach to quantitative
	imaging of complex geological structures. Indeed, wide-aperture acquisition
	geometries offer the opportunity to develop an accurate background
	velocity model for the subsequent waveform tomography. This is critical,
	because the building of the macromodel remains an open question when
	only near-vertical reflection data are considered.},
  file = {:./bibs/shin_2010.pdf:PDF},
  keywords = {Inverse theory
	
	Seismic tomography
	
	Computational seismology
	
	Acoustic properties
	
	finite-difference methods
	
	thrust belt
	
	traveltime and full waveform inversions
	
	wide-aperture seismic data}
}

@ARTICLE{Shipp2002,
  author = {Shipp, Richard M. and Singh, Satish C.},
  title = {Two-dimensional full wavefield inversion of wide-aperture marine
	seismic streamer data},
  journal = {Geophysical Journal International},
  year = {2002},
  volume = {151},
  pages = {325--344},
  number = {2},
  abstract = {Summary A 2-D full wavefield inversion method is presented for the
	processing of wide-aperture data. The diversity of information contained
	within such datasets may be handled in a complete manner by first
	matching the traveltimes of the main events and then progressing
	to waveform fitting of the data through explicit full wavefield modelling.
	Our wavefield inversion scheme is based upon a finite difference
	solution of the 2-D elastic wave equation in the time distance domain.
	The strength of adopting such an approach is the ability to generate
	all possible wave types within a given 2-D model (multiples, converted
	waves, etc.) and thus to simulate and accurately model complex seismic
	wavefields. The aim of the inversion is to find the 2-D P-wave velocity
	model that minimizes the least squared difference between the observed
	and synthetic data across the full range of offsets. Following extensive
	testing on synthetic data, the wavefield inversion scheme has been
	applied to wide-aperture real marine seismic streamer datasets. We
	present results from the synthetic testing and the wavefield inversion
	of wide-aperture real data out to 12 km offset that was recorded
	on a single streamer. Even though current computational restrictions
	allow only a small subsection of the data to be analysed, these examples
	demonstrate the potential value of wide-aperture 2-D full wavefield
	inversion.},
  doi = {10.1046/j.1365-246X.2002.01645.x},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/shipp2002.pdf:PDF},
  issn = {1365-246X},
  keywords = {finite-difference methods, inversion, seismic modelling, seismic structure,
	tomography, waveform analysis},
  owner = {franciane},
  publisher = {Blackwell Science Ltd},
  timestamp = {2011.02.22},
  url = {http://dx.doi.org/10.1046/j.1365-246X.2002.01645.x}
}

@ARTICLE{sirgue2004,
  author = {Sirgue, Laurent and Pratt, R. Gerhard},
  title = {Efficient waveform inversion and imaging: A strategy for selecting
	temporal frequencies},
  journal = {Geophysics},
  year = {2004},
  volume = {69},
  pages = {231-248},
  number = {1},
  abstract = {Prestack migration and/or inversion may be implemented in either the
	time or the frequency domain. In the frequency domain, it is possible
	to discretize the frequencies with a much larger sampling interval
	than that dictated by the sampling theorem and still obtain an imaging
	result that does not suffer from aliasing (wrap around) in the depth
	domain. The selection of input frequencies can be reduced when a
	range of offsets is available; this creates a redundancy of information
	in the wavenumber coverage of the target. In order to optimize the
	use of this information, we define a new discretization strategy
	that depends on the maximum effective offset present in the surface
	seismic survey: the larger the range of offsets, the fewer frequencies
	are required. The strategy, exact in a homogeneous 1D earth, selects
	frequencies by making use of the well-known effect of image stretch
	in normal-moveout (NMO) correction and in migration (usually considered
	detrimental for the imaging). The strategy is also useful in more
	general earth models: we apply it to the 2D Marmousi model and recover
	a continuous range of wavenumbers using only three input frequencies.
	The Marmousi inversion result accurately predicts all other data
	frequencies, demonstrating the redundancy of the data.},
  doi = {10.1190/1.1649391},
  eprint = {http://geophysics.geoscienceworld.org/cgi/reprint/69/1/231.pdf},
  file = {:./bibs/Sirgue_Pratt2004.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/69/1/231}
}

@ARTICLE{stekl1998,
  author = {Stekl, I. and Pratt, R. G.},
  title = {Accurate viscoelastic modeling by frequency-domain finite differences
	using rotated operators},
  journal = {Geophysics},
  year = {1998},
  volume = {63},
  pages = {1779-1794},
  number = {5},
  abstract = {The viscoelastic wave equation is an integro-differential equation
	that requires special methods when using time-domain numerical finite-difference
	methods. In the frequency domain, the integral terms are easily represented
	by complex valued elastic media properties. There are further significant
	advantages to using the frequency domain if the forward of the inverse
	problem requires modeling or inverting a large number of prestack
	source gathers. Numerical modeling is expensive for seismic data
	because of the large number of wavelengths typically separating sources
	from receivers, which results in a need for a large number of grid
	points. A major obstacle to using frequency-domain methods is the
	consequent storage requirements. To reduce these, we maximize the
	accuracy and simultaneously minimize the spatial extent of the numerical
	operators. We achieve this by extending earlier published methods
	introduced for the viscoacoustic case to the viscoelastic case. This
	requires the formulation of two new numerical operators: a differencing
	operator in a rotated coordinate frame and a lumped mass term. The
	new operators are combined with ordinary second-order, finite-difference
	operators in an optimal manner to minimize numerical errors without
	increasing the size of the numerical operator. For a fixed number
	of grid points, the resulting second-order differencing scheme is
	no more expensive than an ordinary second-order differencing scheme,
	but a numerical dispersion analysis shows that the number of grid
	points required per smallest wavelength is reduced from approximately
	15 to approximately 4. The new scheme is also capable of handling
	embedded fluid layers without instability. We demonstrate that no
	further improvement in performance can be achieved using higher order
	spatial operators because of the associated computational overheads
	associated with the larger differencing operators. The new viscoelastic
	modeling scheme is used to study a crosshole data set in which the
	exact nature of the seismic coda is unclear. The results of the modeling
	study indicate this coda is likely related to the generation of mode-converted
	shear waves within the complicated, finely layered sediments at the
	site.},
  doi = {10.1190/1.1444472},
  file = {:./bibs/steklpratt1998.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/63/5/1779}
}

@ARTICLE{Symes2009,
  author = {W W Symes},
  title = {The seismic reflection inverse problem},
  journal = {Inverse Problems},
  year = {2009},
  volume = {25},
  pages = {123008},
  number = {12},
  abstract = {The seismic reflection method seeks to extract maps of the Earth's
	sedimentary crust from transient near-surface recording of echoes,
	stimulated by explosions or other controlled sound sources positioned
	near the surface. Reasonably accurate models of seismic energy propagation
	take the form of hyperbolic systems of partial differential equations,
	in which the coefficients represent the spatial distribution of various
	mechanical characteristics of rock (density, stiffness, etc). Thus
	the fundamental problem of reflection seismology is an inverse problem
	in partial differential equations: to find the coefficients (or at
	least some of their properties) of a linear hyperbolic system, given
	the values of a family of solutions in some part of their domains.
	The exploration geophysics community has developed various methods
	for estimating the Earth's structure from seismic data and is also
	well aware of the inverse point of view. This article reviews mathematical
	developments in this subject over the last 25 years, to show how
	the mathematics has both illuminated innovations of practitioners
	and led to new directions in practice. Two themes naturally emerge:
	the importance of single scattering dominance and compensation for
	spectral incompleteness by spatial redundancy.},
  file = {:./bibs/symes2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.17},
  url = {http://stacks.iop.org/0266-5611/25/i=12/a=123008}
}

@ARTICLE{tabei2002,
  author = {Makoto Tabei and T. Douglas Mast and Robert C. Waag},
  title = {A k-space method for coupled first-order acoustic propagation equations},
  journal = {The Journal of the Acoustical Society of America},
  year = {2002},
  volume = {111},
  pages = {53-63},
  number = {1},
  doi = {10.1121/1.1421344},
  file = {:./bibs/tabei_mast_waag_2001.pdf:PDF},
  keywords = {ultrasonic propagation; differential equations; inhomogeneous media},
  publisher = {ASA},
  url = {http://link.aip.org/link/?JAS/111/53/1}
}

@BOOK{tarantola,
  title = {Inverse problem theory and methods for model parameter estimation},
  publisher = {SIAM},
  year = {2005},
  author = {Albert Tarantola},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Inverse_Problem_Theory_and_Methods_for_Model_Parameter_Estimation.pdf:PDF},
  keywords = {inverse problem},
  owner = {franciane},
  timestamp = {2011.09.05}
}

@ARTICLE{tarantola1984,
  author = {Tarantola, Albert},
  title = {Inversion of seismic reflection data in the acoustic approximation},
  journal = {Geophysics},
  year = {1984},
  volume = {49},
  pages = {1259-1266},
  number = {8},
  abstract = {The nonlinear inverse problem for seismic reflection data is solved
	in the acoustic approximation. The method is based on the generalized
	least-squares criterion, and it can handle errors in the data set
	and a priori information on the model. Multiply reflected energy
	is naturally taken into account, as well as refracted energy or surface
	waves. The inverse problem can be solved using an iterative algorithm
	which gives, at each iteration, updated values of bulk modulus, density,
	and time source function. Each step of the iterative algorithm essentially
	consists of a forward propagation of the actual sources in the current
	model and a forward propagation (backward in time) of the data residuals.
	The correlation at each point of the space of the two fields thus
	obtained yields the corrections of the bulk modulus and density models.
	This shows, in particular, that the general solution of the inverse
	problem can be attained by methods strongly related to the methods
	of migration of unstacked data, and commercially competitive with
	them.},
  doi = {10.1190/1.1441754},
  file = {:./bibs/tarantola1984.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/49/8/1259}
}

@ARTICLE{Valenciano2009,
  author = {Alejandro A. Valenciano and Biondo L. Biondi and Robert G. Clapp},
  title = {Imaging by target-oriented wave-equation inversion},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCA109-WCA120},
  number = {6},
  doi = {10.1190/1.3250267},
  file = {:./bibs/valenciano2009.pdf:PDF},
  keywords = {geophysical signal processing; geophysical techniques; Hessian matrices;
	least squares approximations; wave equations},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/WCA109/1}
}

@ARTICLE{Vigh2009,
  author = {Vigh, Denes and Starr, E. William and Kapoor, Jerry},
  title = {Developing Earth models with full waveform inversion},
  journal = {The Leading Edge},
  year = {2009},
  volume = {28},
  pages = {432-435},
  number = {4},
  abstract = {Exploration in more geologically complex areas requires new methodologies.
	In its quest to answer these new challenges, the oil and gas industry
	has moved from ray-based imaging to finite-difference, wave-equation
	migration to achieve better subsurface descriptions of target zone
	and reservoirs. Notable in this progression is the movement from
	ray-traced Kirchhoff algorithms through one-way wave-equation methods
	to use the acoustic two-way wave equation.},
  doi = {10.1190/1.3112760},
  eprint = {http://tle.geoscienceworld.org/cgi/reprint/28/4/432.pdf},
  file = {:./bibs/vigh2009.pdf:PDF},
  owner = {franciane},
  timestamp = {2011.02.17},
  url = {http://tle.geoscienceworld.org/cgi/content/abstract/28/4/432}
}

@ARTICLE{virieux1986,
  author = {Virieux, Jean},
  title = {P-SV wave propagation in heterogeneous media; velocity-stress finite-difference
	method},
  journal = {Geophysics},
  year = {1986},
  volume = {51},
  pages = {889-901},
  number = {4},
  abstract = {I present a finite-difference method for modeling P-SV wave propagation
	in heterogeneous media. This is an extension of the method I previously
	proposed for modeling SH-wave propagation by using velocity and stress
	in a discrete grid. The two components of the velocity cannot be
	defined at the same node for a complete staggered grid: the stability
	condition and the P-wave phase velocity dispersion curve do not depend
	on the Poisson's ratio, while the S-wave phase velocity dispersion
	curve behavior is rather insensitive to the Poisson's ratio. Therefore,
	the same code used for elastic media can be used for liquid media,
	where S-wave velocity goes to zero, and no special treatment is needed
	for a liquid-solid interface. Typical physical phenomena arising
	with P-SV modeling, such as surface waves, are in agreement with
	analytical results. The weathered-layer and corner-edge models show
	in seismograms the same converted phases obtained by previous authors.
	This method gives stable results for step discontinuities, as shown
	for a liquid layer above an elastic half-space. The head wave preserves
	the correct amplitude. Finally, the corner-edge model illustrates
	a more complex geometry for the liquid-solid interface. As the Poisson's
	ratio v increases from 0.25 to 0.5, the shear converted phases are
	removed from seismograms and from the time section of the wave field.},
  doi = {10.1190/1.1442147},
  file = {:./bibs/virieux1986.pdf:PDF},
  url = {http://geophysics.geoscienceworld.org/cgi/content/abstract/51/4/889}
}

@ARTICLE{Virieux2009,
  author = {J. Virieux and S. Operto},
  title = {An overview of full-waveform inversion in exploration geophysics},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WCC1-WCC26},
  number = {6},
  doi = {10.1190/1.3238367},
  file = {:./bibs/Virieux_2009.pdf:PDF},
  keywords = {geophysical prospecting; seismic waves; seismology; seismometers},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/WCC1/1}
}

@ARTICLE{wang2003,
  author = {Wang, Huazhong and Zhang, Libin and Ma, Zaitian},
  title = {Seismic wave imaging in visco-acoustic media},
  journal = {Science in China Series A: Mathematics},
  year = {2004},
  volume = {47},
  pages = {146-154},
  note = {10.1360/04za0013},
  affiliation = {Tongji University School of Ocean and Earth Science 200092 Shanghai
	China},
  file = {:./bibs/wang_2003.pdf:PDF},
  issn = {1006-9283},
  issue = {0},
  keyword = {Mathematics and Statistics},
  publisher = {Science China Press, co-published with Springer},
  url = {http://dx.doi.org/10.1360/04za0013}
}

@ARTICLE{Wang2009,
  author = {Wang, Yanghua and Rao, Ying},
  title = {Reflection seismic waveform tomography},
  journal = {J. Geophys. Res.},
  year = {2009},
  volume = {114},
  pages = {B03304--},
  number = {B3},
  month = mar,
  abstract = {In seismic waveform tomography, if using reflection data with limited
	source-receiver offsets, it is difficult to reconstruct the deep
	part of the subsurface velocity model. We present two approaches
	to tackle this problem: layer stripping and weighted updating. In
	a layer-stripping procedure, we replace the top portion of seismic
	data with synthetics generated from the previous-layer inversion
	and make the current inversion focus on the minimization of the data
	misfit corresponding to the deep part of the model. To improve efficiency,
	we use only sparsely sampled frequency data in the deeper-layer inversions,
	unlike the first-layer inversion where we use densely sampled frequency
	data as usual. The sparsely sampled frequencies together have the
	full wave number coverage for effective imaging. Combined use of
	dense and sparse sampling in frequency is a compromise between resolution
	and efficiency, as it reduces the number of iterations needed in
	layer-stripping inversion while still producing a good image. In
	the second scheme, we apply depth-dependent weights to model updates
	in order to improve the convergence in an iterative solution. The
	weighting is inversely proportional to the ray density variation
	along the depth and is mathematically equivalent to the application
	of an inverse Hessian matrix which sharpens the gradient vector for
	model updating. For real seismic data, we transfer point source shot
	records to line source records, by partial amplitude compensation
	and phase adjusting, before inputting it to the waveform tomography.
	We perform traveltime inversion to generate a reliable layered velocity
	model and then waveform tomography to produce a high-resolution image
	of the subsurface model through frequency domain iteration.},
  issn = {0148-0227},
  keywords = {exploration seismology, waveform tomography, reflection seismics,
	7270 Seismology: Tomography, 0902 Exploration Geophysics: Computational
	methods: seismic, 0910 Exploration Geophysics: Data processing},
  owner = {franciane},
  publisher = {AGU},
  timestamp = {2011.02.22},
  url = {http://dx.doi.org/10.1029/2008JB005916}
}

@ARTICLE{qingzeng2001,
  author = {Yan Qing Zeng and Qing Huo Liu},
  title = {A staggered-grid finite-difference method with perfectly matched
	layers for poroelastic wave equations},
  journal = {The Journal of the Acoustical Society of America},
  year = {2001},
  volume = {109},
  pages = {2571-2580},
  number = {6},
  doi = {10.1121/1.1369783},
  file = {:./bibs/QingZeng_2001.pdf:PDF},
  keywords = {buried object detection; acoustic applications; finite difference
	methods; acoustic wave propagation},
  publisher = {ASA},
  url = {http://link.aip.org/link/?JAS/109/2571/1}
}

@BOOK{Zhadanov,
  title = {Geophisical Inverse Theory and Regularization Problems},
  publisher = {Elsevier},
  year = {2002},
  author = {M.S. Zhadanov},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/livros/Geophysical_Inverse_Theory_and_Regularization_Problems.pdf:PDF},
  keywords = {inverse problem},
  owner = {franciane},
  timestamp = {2011.09.05}
}

@ARTICLE{Zhou2009,
  author = {Bing Zhou and Stewart Greenhalgh},
  title = {On the computation of the Frèchet derivatives for seismic waveform
	inversion in 3D general anisotropic, heterogeneous media},
  journal = {Geophysics},
  year = {2009},
  volume = {74},
  pages = {WB153-WB163},
  number = {5},
  doi = {10.1190/1.3123766},
  file = {:/home/franciane/Documents/Doutorado/bibliografia/bibs/zhou2009.pdf:PDF},
  keywords = {anisotropic media; inverse problems; Jacobian matrices; perturbation
	theory; seismic waves; seismology},
  owner = {franciane},
  publisher = {SEG},
  timestamp = {2011.02.14},
  url = {http://link.aip.org/link/?GPY/74/WB153/1}
}

